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The (f> 4 scalar field theory in three dimensions, prototype for the study of phase transitions, is 
investigated by means of the Hierarchical Reference Theory (HRT) in its smooth cut-off formulation. 
The critical behavior is described by scaling laws and critical exponents which compare favorably 
with the known values of the Ising universality class. The inverse susceptibility vanishes identically 
inside the coexistence curve, providing a first principle implementation of the Maxwell construction, 
and shows the expected discontinuity across the phase boundary, at variance with the usual sharp 
cut-off implementation of HRT. The correct description of first and second order phase transitions 
within a microscopic, non-perturbative approach is thus achieved in the smooth cut-off HRT. 
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I. INTRODUCTION 

The equilibrium properties of classical fluids have been 
extensively studied in the past decades and very success- 
ful theories of both the gas and the liquid states are avail- 
able by now [3, However, a fully satisfactory descrip- 
tion of the liquid-vapor phase transition within a micro- 
scopic liquid state approach has still to come. Although 
the qualitative features and the universal properties char- 
acterizing the transition are known since the development 
of the Renormalization Group (RG) approach [3[ , a quan- 
titative theory able to predict phase boundaries for spe- 
cific fluid models is not available yet: Most of the theories 
developed in liquid-state physics deal mainly with short- 
range correlations and short-wavelength density fluctua- 
tions. That is why they are not able to reproduce, even 
qualitatively, the phenomenology of the critical region. 
Furthermore, these liquid-state theories give an unsatis- 
factory description of the thermodynamics along the first 
order liquid- vapor transition line. The convexity of the 
free energy is not guaranteed in the integral equation 
approach or in mean-field theories and in fact is always 
violated inside the two-phase region. This intrinsic de- 
ficiency is usually overcome by means of ad hoc proce- 
dures, like Maxwell double tangent construction. An ac- 
curate description of the long-wavelength fluctuations is 
implemented within the liquid-state framework in the Hi- 
erarchical Reference Theory of fluids (HRT) [4|, |5[ , where 
the basic RG concept of selective turning on of fluctu- 
ations is built into a genuine liquid-state theory. This 
feature made HRT the only practical scheme able to de- 
termine both universal and non-universal critical proper- 
ties of fluids, so far. Moreover, the proper treatment of 
long- wavelength fluctuations forces the free energy to dis- 
play the correct convexity also in the two-phase region, 
thereby incorporating Maxwell construction: A major 
achievement of HRT. Although quite successful, the cur- 
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rent sharp cut-off formulation of HRT predicts a diverg- 
ing compressibility along the binodal (in any dimension 
smaller than four) Q while, for a scalar order parameter, 
a finite value of the inverse compressibility at coexistence 
is expected. In this paper we show that this failure of 
the HRT approach in the sharp cut-off formulation can 
be overcome by studying the smooth cut-off implemen- 
tation, already introduced several years ago [7| . In order 
to clarify the qualitative properties expected within this 
approach, we consider the simplest model displaying a 
phase transition within the Ising universality class: The 
scalar </> field theory in three dimensions. A somehow re- 
lated formalism has been recently put forward by Caillol 
_$] who developed a non-perturbative renormalization- 
group theory for fluids, within the Grand Canonical en- 
semble and applied it to the Kac model. 

The paper is organized as follows: In Section II we 
review the sharp and the smooth cut-off formalism for 
a general microscopic fluid model. In Section III this 
formalism is applied to a 4 field theory, where several 
simplifications occur. An approximate closure relation 
is introduced and discussed. The numerical solution of 
the equations at low temperatures and the comparison 
between the two approaches is presented in Section IV, 
where the behavior at coexistence is analyzed. The crit- 
ical properties of HRT in the smooth cut-off formulation 
are derived in Section V, while Section VI contains a 
short summary of the results and some perspective. 

II. HRT EQUATIONS IN SHARP AND 
SMOOTH CUT-OFF FORMULATIONS 

The common principle underlying HRT, in its various 
formulations, is the selective treatment of density fluctu- 
ations on different length scales. Starting from a mean 
field approximation, the effects of density fluctuations 
on short wavelengths are included first. This is accom- 
plished by turning on the attractive part of the potential 
progressively, starting from its Fourier components with 
larger wave vectors. A key role is therefore played by the 
cut-off wave vector Q which separates the Fourier com- 
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ponents already taken into account (those with k > Q) 
from those ignored (k < Q). The full spectrum of fluctu- 
ations is then included only in the Q — > limit, while at 
Q — > oo the mean field approximation is recovered. The 
variation in the free energy induced by an infinitesimal 
change in Q can be computed exactly by use of perturba- 
tion theory, leading to a differential "evolution" equation 
for the thermodynamics of the system. 

We study a system of particles in dimension d, inter- 
acting via a two-body potential v(r) formally written as 
the sum of two contributions: 



v(r) = vji{r) + w(r) 



(1) 



where Vn(r) is a short-range, repulsive reference part, 
and w(r) is a mostly attractive term defined by a regular 
function which can be Fourier transformed. The proper- 
ties of the system interacting via vr are assumed to be 
known by other liquid state theories or simulations. In or- 
der to implement the principle of HRT, the Fourier com- 
ponents of w(r) have to be included gradually. Therefore 
we define a sequence of Q systems characterized by an in- 
teraction vq(t): 



VQ (r) = v R (r) + w Q (r) 



(2) 



interpolating between the reference system (obtained for 
wq (r) = 0) and the fully interacting model (for wq (r) — 
w(r)). The difference between sharp and smooth cut-off 
formulation resides in the way the cut off on the Fourier 
components of the potential is implemented, as discussed 
in the following. 



A. Sharp cut-off 



vector since the very beginning of the integration proce- 
dure. The k = term in fact gives rise to the mean-field 
contribution to the free energy which provides a physi- 
cally correct starting point for the description of phase 
transitions. This is achieved by introducing a modified 
free energy density defined as: 



A 
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= 0) - <fo(fc = 0) (4) 



where </>(r) = —j3w{r), [3 — and p is the par- 

ticle density. Analogously, the modified direct correla- 
tion function includes the Random Phase Approximation 
(RPA) contribution: 



C Q (k) = CQ(k)+<f>(k)-<fi Q (k) 



(5) 



where cq (r) is the usual direct correlation function [l| of 
the Q-system (with the inclusion of the ideal gas term). 
The evolution of the (modified) free energy satisfies the 
exact differential equation: 



dA 
dQ 
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Cq(Q) 



(G) 



where Qd is the volume of the unit sphere in d dimensions 
divided by (27r) d . An analogous equation governs the 
evolution of Cq(/c) which, however, depends on the three 
and four body correlation functions of the Q-system, 
thereby giving rise to an exact hierarchy of differential 
equations. In the following Section we will provide an 
approximate closure at the level of the first equation of 
the hierarchy, Eq. ©. 



In the sharp cut-off formulation, the precise definition 
of wq is conveniently given in Fourier space: 



w Q {k) 



w(k) for k > Q 
for k < Q 



(3) 



where the tilde denotes Fourier transform. Here, the cut- 
off Q varies between Q — oo, where the potential reduces 
to the reference part, and Q = 0, where the physical 
interaction (fT]) is recovered. Due to the long range os- 
cillating tail in u>Q(r), induced by the sharp cut-off, no 
phase transition occurs in each Q-system until the Q — > 
limit is attained. It is now possible to write a differential 
equation for the evolution of the free energy Aq of the Q- 
system when the cut-off is varied. This equation, though 
exact, is not closed because it involves the two body cor- 
relation function of the Q-system, and suitable approx- 
imations linking thermodynamics and correlations must 
be introduced. Before discussing the adopted closure, 
following the RG approach, it is convenient to introduce 
an additional definition in order to take into account the 
Fourier component of the potential with vanishing wave 



B. Smooth cut-off 



In the smooth cut-off formulation of HRT Q the dis- 
continuity in the Fourier transform of wq induced by the 
definition ([3]) is removed by the alternative choice 



w t (k)=w(k)~e 2t w{ke t ) 



(7) 



which identically vanishes at t = and tends to w(k) 
for t — > oo. If w(k) varies on the characteristic wave 
vector scale A, the quantity Q — Ae~* acts as an ef- 
fective cut-off wave vector: The Fourier components of 
wt{k) for k < Q are efficiently suppressed and the long 
range repulsive tail present in w t (r) contrasts the ten- 
dency towards phase separation at any finite t. The se- 
quence of intermediate potentials wAk) defined in (JT]) 
belongs to the class studied in Ref. [8| and corresponds 
to the particular choice of cut-off function Rt{k) w(k) = 



e- 2 *w}(fce*i[M)(fc) 



-2t„r, 



(fee 4 )]" 



The requirements 



stated in [8| are satisfied by Rt(k) provided w(k) decays 
rapidly at large wave vectors. First order perturbation 
theory provides the change in the free energy induced by 
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an infinitesimal increase of the parameter t. By use of 
the same the definition (0J already introduced, the exact 
evolution equation of the (modified) free energy in the 
smooth cut-off formulation of HRT reads: 



dA t 
dt 



d d k p( ,,d4>t(k) 



(2tt) 



dt 



(8) 



where </>t(fc) = —(3wt(k) and the two point function 
F t (k) = — [c f (k)]~ is just p times the structure factor 
Sq(k) of the Q-system. As usual, it is more convenient 
to introduce the modified direct correlation function Ct 
(Eq. ([5])) which includes the total interparticle poten- 
tial in mean-field approximation, because we expect that 
its convergence properties as t — > oo will be better than 
those of the "bare" correlations c t (k). By use of Eqs. ([5]) 
and |(7j) we get a formal expression for the bare two point 
function appearing in Eq. (jSJ: 



Ft{k) 



Ct(k) - 4>{k) + fa(k) 
C t (k) -e~ 2t 0(fce*) 



(9) 



Analogously to the sharp cut-off case, the flow equation 
(j8]) for At is not closed and a suitable approximation for 
Ct(k) must be introduced. 



III. HRT EQUATIONS FOR THE 
FIELD THEORY 



SCALAR 



In order to test the accuracy of HRT in the description 
of phase transitions, we specialize the general equations 
<j6j) and ijHJ) to a particular model of fluid. We first note 
that, in systems with regular short-range attractive inter- 
actions, the Fourier transform of 4>{r) is a positive, mono- 
tonic decreasing function of the wave vector, quadratic in 
k at small k. The simplest model with these features is 
defined by an effective interaction of the form: 



b (A 2 




k 2 ) for k < A 
for k > A 



(10) 



By substituting this form into Eqs. 
respectively 



© and © we get, 



dA 



dQ 
dAq 
dQ 



= bQ 
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k<Q (2^C Q (fc)+6(fc 2 -Q 2 ) 



(12) 



where the cut-off wave vector in the smooth cut-off equa- 
tion (fT2"|) is defined as Q = Ae~*. The equations hold for 
Q < A and must be supplemented by the initial con- 
dition at Q — A, where Aq reduces to the mean field 
expression. In both cases, as initial condition, we take 



a Landau- Ginzburg quartic form, centered around the 
critical density p c , with coupling constants (r, u) [3[: 

A A (p) = A A (Pc) ~r{p- Pc) 2 -u(p- pc) 4 (13) 

In a fully microscopic model, defined by a realistic two 
body interaction w(r), we expect that the quartic form 
(|13[) represents the physical free energy density in a small 
neighborhood of the critical density p c . In this case, the 
effective couplings (r, u) are functions of the tempera- 
ture and may be estimated in mean field theory. Short 
range density fluctuations, corresponding to wave vec- 
tors Q > A, are expected to renormalize the bare values, 
as shown in Ref. Q where this problem is examined 
within the HRT formalism. The modified direct cor- 
relation function appearing in the two evolution equa- 
tions (fTTjl . (fT2")l is approximated in the spirit of the Local 
Potential Approximation of RG [1, 0], by an Ornstein- 
Zernike form 



Co(fc) = C Q (0) - bk 2 = 



d 2 A Q 
dp 2 



bk 2 



(14) 



In Eq. (|14p the range of the direct correlation function 
has been assumed to coincide with that of the chosen 
potential (|10| with no renormalization due to fluctua- 
tions. This is actually the only approximation introduced 
into the HRT formalism and implies an analytic momen- 
tum dependence of the correlation functions in the whole 
phase diagram, critical point included. In order to allow 
for a non vanishing critical exponent r\ within HRT, we 
must go beyond the parametrization (|14[) taking into ac- 
count the effects of fluctuations on the range of the corre- 
lation function. This can be accomplished by examining 
the second equation of the HRT hierarchy, governing the 
evolution of Cq(A;) Q]. In fact, it has been shown 0] that, 
in the critical region and in the limit Q — > 0, this equation 
allows one to reproduce correctly the critical exponents 
up to the e 2 term in the e = 4 — d expansion, hence giv- 
ing a non- vanishing value of rj. The last equality in (|14|) 
directly follows from the compressibility sum rule [l[ , ac- 
cording to which the structure factor evaluated at zero 
wave vector is equal to the reduced compressibility of the 
system. By use of this parametrization, Eqs. (fTTj) . (fT2]) 
give rise to a couple of partial differential equations de- 
scribing the effect of fluctuations into the mean field free 
energy within the sharp and smooth cut-off formulation 
of HRT for a potential of the form (jTOjl . The resulting 
HRT equations are then: 



dAq 
dQ 

dA Q 
dQ 



rf-i 



In 



-A'/ 



bA 2 



-Aq + bQ 2 



bil d Q 



d+l 



-A% + bQ 2 



(15) 
(16) 



for the sharp and smooth cut-off respectively, where 
primes mean differentiation with respect to p. The sharp 
cut-off equation (fTS"]) for a (/> 4 field theory has already 
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been studied in Refs. [H, Q both in the critical region 
and at phase coexistence where, for Q — » 0, the term 
— Aq + bA 2 in the argument of the logarithm can be 
neglected. We remark that this does not affect by any 
means the universal behavior at criticality or the qualita- 
tive features of the first-order transition discussed here. 
In Ref. [l(| the smooth cut-off RG equations have been 
derived and studied for a field theory. The coinci- 
dence, after a trivial rescaling, of our HRT equation (fTB]) 
with Eq. (2.6) of Ref. [l(| for the special cut-off choice 
n = | + 1 shows that indeed i) the parabolic potential 
model defined by (fTU)) is equivalent to a d 4 scalar field 
theory and ii) that the smooth cut-off HRT formalism 
becomes equivalent to the RG in the scaling limit [ill ]. 



IV. THE COEXISTENCE BOUNDARY 

A fully implicit predictor-corrector finite difference al- 
gorithm is used in order to solve the HRT evolution equa- 
tions (|I5II6[) for a <f> A field theory. The initial condition 
(fT5| is imposed at Q — A, consistently with the assumed 
form of the interaction potential (fit))) . The inverse com- 
pressibility x 1 — —d 2 Ag/dp 2 at the end of integration 
(i.e. Q — ► 0) is shown in Fig. [T] as a function of p — p c 
for a representative choice of parameters in the broken- 
symmetry regime. At large densities the two formula- 
tions of HRT provide very similar results. Moreover, both 
show a region of infinite compressibility, consistent with 
the convexity requirement of the free energy. The impor- 
tant difference between the smooth and the sharp cut-off 
formulation is the presence of a discontinuity across the 
coexistence curve in three dimensions in the smooth cut- 
off case, in agreement with the expected behavior for a 
scalar order parameter. Instead, as already discussed in 
Ref. 6], the sharp cut-off HRT equation predicts the 
divergence of the compressibility when coexistence is ap- 
proached. 

We now provide an analytical interpretation of the ori- 
gin of the discontinuity following the procedure devised 
in [||, for the sharp cut-off formalism. We first write our 
equation (|1€5[) in a scale invariant form. By applying the 
rescaling: 




Aq 



to Eq. (TT6|) . we obtain: 

~ df Q d+2 



dQ J xx + Q 2 



(17) 



(18) 



d 2 f 

where we have set f xx = -g^i- By performing the substi- 
tution 



1 



fxx H~ 



(19) 




FIG. 1: Inverse compressibility as a function of the order 
parameter p — p c in three dimensions from the sharp (above) 
and smooth cut-off formulation (below). Units are set by 
6A 2 = 1 and 5^3 A 3 = 1, while the coupling constants are: 
r = —0.47 and u = 0.0035. Upper panel: result from the 
sharp cut-off formulation, Eq. (|15[) . Lower panel: result from 
the smooth cut-off formulation, Eq. (|16[1 . Notice the presence 
of the expected discontinuity across the coexistence boundary 
in the smooth cut-off case. 



and taking the double derivative of Eq. (|18[) with respect 
to x, the equation is written in quasi-linear form, suitable 
for an analytical study: 



1 86 



-Q 



(20) 



In the two-phase region, the observed convexity of the 
free energy implies the divergence of the compressibility 
X- This means that for Q — > 0, f xx °c X -1 vanishes and 
6 — > +oo. In this case, we can neglect the first term of 
Eq. ([20]) . so that: 



= -2Q 



whose solution is 



x 2 )Q~ 



(21) 



(22) 



where the integration constant a^o plays the role of 
rescaled coexistence density. Indeed, this expression 
shows that 6 — > +oo as Q 0, consistently with our 
assumption, only for \x\ < xq: HRT correctly predicts 
the existence of a finite region of infinite compressibility, 
but our analysis is not able to describe the behavior of 6 
across the phase boundary, i.e. the transition between a 
finite solution for > xq and the asymptotic form (|2"2")) 
inside the binodal. Following Q, it is useful to zoom-in 
the region close to Xq by rescaling the x variable as: 



z — (x — x + aQ 2 )Q 



(23) 
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The additive term aQ 2 takes into account the fluctuation 
corrections to the position of the phase boundary, which 
are known to play an important role in determining the 
asymptotic solution on the binodal [||. Then Eq. (j2"0|) 
becomes: 



1 



oQ az 



-Q 



+ 2Q 2 

dz 2 



(24) 



By keeping only the dominant terms as Q — > in Eq. 
I, we obtain the following fixed-point equation: 



(25) 



2adj)_ _ d 2 (j) 
(f> 2 dz dz 2 

Once integrated, the equation gives: 

<p -\ In (0c — 2a) = — c; 

c 



(26) 



where c, k are integration constants. Inside the coexis- 
tence region, i.e. for z — > — oo, Eq. (|26[) must match Eq. 
which, expressed in terms of z, reads 



<f> 



-(at + xq)z w — 2xoz 



(27) 



This sets the constant c at the value c = 2xo . When 
the coexistence boundary is approached from outside, i.e. 
for z — > +oo, Eq. (|26[) predicts a finite compressibility: 
<f> — > 2a/c = 0o- We can then rewrite Eq. (|26[) as: 



<fi + 4>o ln((/>/</>o — 1) = — 2.tq^ + const 



(28) 



which describes the asymptotic behavior of the solution 
across the phase boundary. A numerical verification of 
such a scaling law can be obtained by considering the 
derivative of Eq. (j28|) with respect to z: 



1 



dz 2xq <j> — 



1 



(29) 



The validity of this equation has been tested on the nu- 
merical solution of Eq. I|16p by selecting small Q values 
(Q < 2 x 10~ 2 ) and mesh points extremely close to the 
phase boundary, so that the rescaled variable z given 
in Eq. (j2"5|) is 0(1). The results are shown in Fig. [21 
where several sets of data points, corresponding to dif- 
ferent Q's and different x's, are shown. The asymptotic 
collapse of the data towards unity, as required by Eq. 
(|29p. can be appreciated. The results shown in the fig- 
ure were obtained with a mesh Ap — 4 x 10 -7 . Probing 
the asymptotic regime further would have required an 
even smaller Ap, so as to sample the region z = 0(1) 
for Q smaller than those considered in the figure. The 
discontinuity of across the binodal, shown in Fig. 1, 
is therefore a genuine consequence of the smooth cut-off 
HRT formalism. The analysis has been performed for a 
particular choice of the potential (JTTJJ) and for a closure 
of the form (fl"4"]) . because these assumptions considerably 
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FIG. 2: Numerical results for the quantity * of Eq. [29] Dif- 
ferent lines correspond to different Q values. From top to 
bottom: InQ = —4.0 (solid line), InQ = —4.2 (dashed line), 
InQ = -4.4 (dotted line), InQ = -4.6 (circles), InQ = -4.8 
(squares), InQ = —5.0 (triangles), InQ = —5.2 (stars), 
InQ = —5.4 (diamonds), InQ = —5.6 (crosses). 



simplify the evolution equation (fTS)) . We also studied the 
smooth cut-off HRT equation for a realistic Yukawa inter- 
action w(r), allowing for a difference between the range 
of the potential and that of the direct correlation func- 
tion. The numerical solution in the asymptotic Q — » 
regime shows the expected jump in the inverse suscepti- 
bility across the coexistence boundary. Our findings are 
consistent with the numerical analysis of Ref . [l(| , where 
the behavior of x has been investigated for a one- 
parameter family of smooth cut-off RG theories which 
includes, as a limiting case, also the sharp cut-off form 
(fl"5)) . besides our Eq. (To) : the discontinuity is present 
for all values of the parameter and vanishes only in the 
sharp cut-off limit. 



THE CRITICAL REGION 



In Section III we mentioned that the smooth cut-off 
HRT equation for a parabolic potential ((TO)) reduces to 
a known RG structure, analogously to the sharp cut-off 
case [ij . It is therefore interesting to evaluate the critical 
exponents, by studying the fixed point solution of the 
HRT equation and compare with the direct numerical 
solution of Eq. fTo]) . Following the usual HRT procedure 
0] , inspired by the RG approach, we first rescale the 
variables At(p) and p (at fixed temperature T), in order 
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to blow-up the critical region: 

H(y,Q) = ^Q- d [A Q { Pc ) - A Q {p)\ 



V 



b 

crQ 2 (P-Pc) 



(30) 



where p c is the critical density. By performing the sub- 
stitution we get: 



Q 



dH d-2 



dQ 



yH' - dH 



1 



1 



1 W 



1 



(31) 



where Hq is the the second derivative of H(y) evaluated 
at the critical density. We first look for a fixed point so- 
lution of Eq. (|3ip i.e. a solution H*(y) independent of 
Q. According to our definition, H(0, Q) — and then 
also H*(Q) = 0. The fixed point solution is an even func- 
tion of y, so that H* (0) = while H* (which explicitly 
appears in Eq. (|3Tj)) has to be tuned in such a way that 
the fixed point solution H*(y) is regular on the whole 
real axis. Numerically it is more convenient to write the 
equation for the quantity = H* (y), which can be 
obtained by differentiation of Eq. (|3"Tj) with respect to y: 



d-2 



-yp 



d + 2 



■ P = 



P 



{p 1 + I) 2 



(32) 



with initial conditions fi(0) = and /i'(0) = p. Equa- 
tion ([32]) admits a single regular solution for a unique 
choice of p (besides the trivial gaussian solution (J,(y) = 
which corresponds to p = 0). By numerical integration 
of the fixed point equation (|32p using a trial and error 
method we found a regular solution for p = —0.186 .... 
By setting H(y,Q) = H*{y) 4- h(y)Q~ x and linearizing 
the evolution equation (f3"Tj) for small h, we obtain the 
eigenvalue equation: 



d-2 



yti + (\-d)h 



h" 



h" 



(33) 



where again h!^ is the second derivative of the eigenfunc- 
tion evaluated at y = 0. By shifting the eigenfunction 
h(y) — > h(y) + const, the constant term in Eq. (|33|) can 
be eliminated. The shifted function satisfies the bound- 
ary conditions h(0) = 1 and h'(Q) = 0. The eigenvalue A 
is determined by the requirement that h(y) is free from 
singularities at all y's. Relevant perturbations, driving 
the system out of the critical point, correspond to pos- 
itive eigenvalues. All the eigenfunctions can be classi- 
fied according to their parity with respect to the sym- 
metry y — > —y: Even solutions describe perturbations 
along the critical isochores, i.e. governing the tempera- 
ture dependence of the free energy in the critical region, 
while odd eigenfunctions correspond to changes along the 
density axis. As usual, a single relevant odd eigenfunc- 
tion is present: h(y) = H* (y) with associated eigen- 
value A = (d - 2)/2 H which implies, via scaling rela- 
tions, a critical exponent r\ = 0, in agreement with the 
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FIG. 3: A log- log plot of the inverse dimensionless compress- 
ibility as obtained by numerical integration of Eq. (|16[l as 
a function of the reduced "temperature" t = \r — r c \. Tri- 
angles refer to the critical isochore, above the critical tem- 
perature r > r c . Circles correspond to data along the co- 
existence curve for r < r c . A linear fit gives an exponent 
7 = 1.28 ± 0.02 close to the result obtained from the fixed 
point analysis 7 = 1.300 ± 0.001 (solid line). 



assumed analyticity of the correlation function (|14|) . The 
numerical solution of the eigenvalue equation (f33|) shows 
that a single relevant even eigenfunction is present, cor- 
responding to an eigenvalue A = 1.539 ± 0.001 in three 
dimensions. The exponent 7 governing the divergence 
of the compressibility x ~ |T — T c |~ 7 is then given by 
7 = 2/A = 1.300 ± 0.001, to be compared with the ac- 
cepted value for the Ising universality class 7 ~ 1.24. 
This result agrees, within numerical uncertainties, with 
the estimate of 7 extracted from the numerical solution 
of the smooth cut-off HRT equation (fTS)) , as shown in 
Fig. [3] Since the usual scaling relations among criti- 
cal exponents are satisfied by HRT and the approximate 
Ornstein-Zernike closure |T5|) implies that the scaling ex- 
ponent 77 is vanishing, knowledge of 7 is sufficient to de- 
termine all the remaining critical exponents. In partic- 
ular, the critical exponent (3 describing the shape of the 
coexistence curve (\p — pA ~ \T — T C \P) is obtained by use 
of the scaling relations |4j, which give (3 = 7/4 ~ 0.325. 
This value is consistent with the fit of the numerical re- 
sults, as shown in Fig. |4j Contrary to the prediction of 
HRT in the sharp cut-off formulation, the smooth cut- 
off equations lead to a divergence of the specific heat 
at criticality: Cy ~ \T — T c \~ a . The critical exponent 
can be evaluated by use of scaling laws with the result 
a = (4 - 37)/2 w 0.050. In Tab. I the values of the 
critical exponents in a three-dimensional fluid according 
to different approximations are reported. "Exact" values 
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FIG. 4: Location of the coexistence curve as obtained by 
numerical integration of Eq. (|16p as a function of the reduced 
"temperature" t = \r — r c \. A linear fit gives f3 = 0.33 ± 0.03, 
consistent with the result from the fixed point analysis f3 ~ 
0.325 (solid line). 



Exponent 


"Exact" 


Sharp 


Smooth 




value 


cut-off HRT 


cut-off HRT 


a 


0.110 


-0.07 


0.05 


P 


0.327 


0.345 


0.330 


1 


1.237 


1.378 


1.300 


S 


4.789 


5 


5 


V 


0.036 









TABLE I: HRT estimates of the critical exponents for the 
Ising universality class in three dimensions compared to 
the exact values [l2l ] obtained by extrapolation of high- 
temperature series expansions. Sharp cut-off HRT results 
from Ref. Q 

are derived from extrapolation of high-temperature se- 
ries expansions [l^] , while the sharp cut-off HRT results 
have been obtained in [f| . In both HRT formulations the 
exponent r\ vanishes because of the OZ assumption, and 
5 = 5 follows from scaling. 

VI. CONCLUSIONS 

The smooth cut-off HRT equation §5§ has been de- 
rived in the framework of liquid state theory for a general 



fluid model. The theory has been then specialized to a 
"parabolic potential" which corresponds to a microscopic 
realization of a <fi 4 theory. The resulting partial differen- 
tial equation turns out to be equivalent to a particular 
formulation of smooth cut-off RG theory [lOj . We stress 
that, contrary to other RG methods, the gradual turning 
on of fluctuations which characterizes the HRT approach, 
is performed directly on the physical quantities, like the 
free energy and no a priori mapping to effective models 
is required. This allows to preserve the information on 
both the universal and the non-universal properties of 
the system. The equation has been numerically solved 
both above and below the critical temperature showing 
that smooth cut-off HRT provides an extremely promis- 
ing tool for the description of the phase diagram of fluids. 
Among the unique features of this formalism we recall 

• A treatment of the critical region consistent with 
scaling laws and characterized by non classical crit- 
ical exponents in good agreement with the exact 
values. 

• A built-in mechanism leading to the convexity of 
the free energy which implements Maxwell's con- 
struction via the inclusion of long range density 
fluctuations. 

• An accurate description of the first order transi- 
tion, which predicts the correct jump of the inverse 
compressibility at the phase boundary. 

Although in this paper we used the smooth cut-off HRT 
for the analysis of a coarse-grained hamiltonian, the <p 
field theory, Eq. ([5]) can be directly applied to fully mi- 
croscopic models of fluids, as already shown in the sharp 
cut-off case (see [4J and references therein). A particu- 
larly favorable choice may be the Yukawa fluid, which 
allows for an exact implementation of the core condition 
for the whole sequence of intermediate models ([7]) inter- 
polating between the reference and the fully interacting 
systems. This feature can make the smooth cut-off for- 
mulation a valuable tool for the theoretical investigation 
of the phase behavior of simple and complex fluids. 
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